to appear on the Proceed, of the Conf. "Beyond the Desert 03" 

Rindberg Castle, Germany, June 2003 

Results from DAMA/Nal and perspectives for 
DAMA/LIBRA 

R. Bernabei\ P. Belli\ F. Cappella\ F. Montecchia^ *, F. Nozzoli\ 

A. Incicchitti^, D. Prosperi^, R. CeruUi^, C.J. Dai*, H.H. Kuang*, J.M. Ma*, 

and Z.P. Ye* ** 

^ Dipartimento di Fisica, Universita di Roma "Tor Vergata" and INFN, Sezione di 

Roma2, 1-00133 Rome, Raly 
^ Dipartimento di Fisica, Universita di Roma "La Sapienza" and INFN, Sezione di 

Roma, F00185 Rome, Raly 
^ INFN - Laboratori Nazionali del Gran Sasso, 1-67010 Assergi (Aq), Italy 
* IHEP, Chinese Academy, P.O. Box 918/3, Beijing 100039, China 

Abstract. The ~ 100 kg highly radiopure Nal(Tl) set-up of the DAMA project 
(DAMA/Nal) took data over seven annual cycles up to July 2002 and has achieved 
results on various rare processes. Its main aim has actually been the exploitation of 
the model independent WIMP annual modulation signature. After this conference the 
total exposure, collected during the seven annual cycles, was released. This cumulative 
exposure (107731 kg x day) has given a model independent evidence for the presence 
of a Dark Matter particle component in the galactic halo at 6.3 a C.L.; this main result 
is summarised here. Some of the many possible corollary model dependent quests for 
the candidate particle are mentioned. At present, after about five years of new devel- 
opments, a second generation low background set-up (DAMA/LIBRA with a mass of 
~ 250 kg Nal(Tl)) was built and is taking data since March 2003. New R&D efforts to- 
ward a possible Nal(Tl) ton set-up, we proposed in 1996, have been funded and started 
in 2003. 



1 Introduction 

The DAMA project has been proposed by the Italian group to INFN and firstly 
funded in 1990 [1]; in 1992 the Chinese colleagues joined the project. DAMA is an 
observatory for rare processes based on the development and use of various kinds 
of radiopure scintillators. Several low background set-ups have been realised; the 
main ones are: i) DAMA/Nal (~ 100 kg of radiopure Nal(Tl)), which took data 
underground over seven annual cycles and was put out of operation in July 2002; 
ii) DAMA/LXe (~ 6.5 kg liquid Xenon), which has been upgraded various times 
and is in operation; iii) DAM A/R&D, which is a set-up devoted to measure- 
ments on prototypes/small-scale-experiments and is in operation; iv) the new 
second generation set-up DAMA/LIBRA (~ 250 kg; more radiopure Nal(Tl)) 
in operation since March 2003. Moreover, in the framework of devoted R&D 
for higher radiopure detectors and PMTs, sample measurements are regularly 

* also: Universita "Campus Bio-medico" di Roma, 00155, Rome, Italy 
** also: University of Zhao Qing, Guang Dong, China 
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carried out by means of the low background DAMA/Ge detector, installed deep 

underground since > 10 years and, in sonic cases, by means of Ispra facilities. 
Recent results from DAMA/LXe are presented in these proceedings elsewhere 
[2], while in ref. [3] some recent results achieved with DAMA/R&D can be found. 

The DAMA/Nal set-up and its performances have been described in details 
in ref. [4]. Since then some upgrading has been carried out; in particular, in 
summer 2000 the electronic chain and data acquisition system were completely 
substituted as reported in ref. [5]. 

As mentioned, DAMA/Nal has taken data over seven annual cycles up to 
July 2002 [5]. The first data release and publication of DAMA/Nal data oc- 
curred in 1996 [6] for an exposure of 4123 kg x day (DAMA/NaI-0) analysed 
by using the pulse shape discrimination (PSD) technique in Nal(Tl) in order 
to discriminate between electromagnetic background and recoils^. Investigation 
of possible diurnal effects has been carried out as well [8]. The main aim of 
the experiment was actually the investigation of the presence of a Dark Mat- 
ter particle component in the galactic halo by means of the model independent 
WIMP annual modulation signature, which has been deeply investigated over 
seven annual cycles (107731 kg x day total exposure). 

DAMA/Nal has also obtained results on the investigation of some exotic 
Dark Matter candidates (such as neutral SIMPs, neutral nuclearities and Q- 
balls) [9], has searched for solar axions [10] and has given results on several 
other rare processes such as e.g. spontaneous emission of protons in ^^Na and 
^^^I with violation of the Pauli exclusion principle [11], nuclear level excitation 
of ^^Na and ^^'^l during charge- non-conserving processes [12], electron stability 
and non-paulian transitions in Iodine atoms by L-shcll[13]. 

In the following; the result obtained by DAMA/Nal in the investigation of the 
WIMP annual modulation signature over the seven annual cycles will be briefly 
summarised. This cumulative 107731 kg x day exposure has given a model in- 
dependent evidence for the presence of a Dark Matter particle component in the 
galactic halo at 6.3 a C.L.; this main result and corollary interpretations in some 
of the possible different model scenarios have been discussed in experimental and 
theoretical details in ref. [5] and have firstly been presented after this Conference 
in a seminar at Gran Sasso National Laboratory (LNGS) at end of July 2003. 

The new second generation larger mass and higher radiopure Nal(Tl) set-up, 
named DAMA/LIBRA, is now in operation and will be briefly introduced. 



^ The Nal(Tl) Dark Matter experiments have discrimination capability as largely 
stated in literature by several groups and at present exploited as main goal by the 
NAIAD experiment. In Zeplin I the pulse shape discrimination technique has been 
employed in LXc as well (sometimes this rejection procedure has boon indicated as 
rejection by timing), similarly as done and published by DAMA/LXe in 1998. More- 
over, we remind that the WIMP annual modulation signature acts (as pointed out 
since the 80's by [7]) itself as a powerful background rejection procedure and that it 
can be effectively exploited only in large mass experiments such as so fax DAMA/Nal 
and DAMA/LIBRA. 
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2 Some generalities 

The DAMA/Nal experiment has been built in order to have suitable mass, sen- 
sitivity and control of the running conditions to investigate the WIMP presence 

in the galactic halo by the annual modulation signature. We remind that the 
WIMP annual modulation signature is based on the annual modulation of the 
signal rate induced by the Earth revolution around the Sun; as a consequence, 
the Earth will be crossed by a larger WIMP flux roughly in June (when its 
rotational velocity is summed to the one of the solar system with respect to 
the Galaxy) and by a smaller one roughly in December (when the two veloci- 
ties arc subtracted) . The annual modulation signature is very distinctive since a 
WIMP-induced seasonal effect must simultaneously satisfy all the following re- 
quirements: the rate must contain a component modulated according to a cosine 
function (1) with one year period, T, (2) and a phase, to, that peaks around ~ 2"'' 
June (3); this modulation must only be found in a well-defined low energy range, 
where WIMP induced recoils can be present (4); it must apply to those events 
in which just one detector of many actually "fires" {single-hit events), since the 
WIMP multi-scattering probability is negligible (5); the modulation amplitude 
in the region of maximal sensitivity is expected to be <7% (6). For the sake 
of completeness, we mention that this latter rough limit would be larger cither 
in case the WIMPs would match the scenario of rcf. [14] (because of kinematic 
effects) or the scenario of ref. [15] (because of a possible external contribution 
to the dark halo from the Sagittarius Dwarf Tidal Stream). To mimic such a 
signature spurious effects or side reactions should not only be able to account 
for the whole observed modulation amplitude, but also to contemporaneously 
satisfy all the requirements; no one has been found or suggested by anyone over 
about a decade. 

As mentioned, to point out the modulation component of a signal, large 
mass apparata with suitable performances and control of the operating condi- 
tions are necessary, such as the ~ 100 kg highly radiopure Nal(Tl) DAMA/Nal 
set-up - which has been the only experiment able to effectively exploit such a 
WIMP signature for about a decade - and now the ~ 250 kg higher radiopure 
Nal(Tl) DAMA/LIBRA set-up. This approach allows to investigate the pres- 
ence of a Dark Matter particle component in the galactic halo independently on 
astrophysical, nuclear and particle physics modelling. Corollary quests for the 
candidate particle require instead to choose a model; therefore, the results are 
not general and refer case by c;ase to the considered model framework ^ as it is 
always the case for results given in form of exclusion plots (which, thus, have 
no generality) and for the specific parameters of a WIMP candidate (mass and 
cross sections) derived from the indirect searches (for some discussion see [5]). 

^ Wc remark that a model framework is identified by the general astrophysical, nuclear 
and particle physics assumptions and by the set of values used for all the experimental 
and theoretical parameters needed in the calculations (for example WIMP nature, 
WIMP couplings, form factors, spin factors, scaling laws, quenching factors, halo 
model, WIMP local velocity, etc., which are affected by relevant uncertainties; see 
for example discussions in ref. [5]). 
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The presence of a model independent positive evidence in the data of DA- 

MA/Nal has been firstly reported by the DAMA collaboration at the TAUP 
conference in 1997 [16] and published also in [17], confirmed in [18,19], further 
confirmed in [20-24] and conclusively confirmed, at end of experiment, in 2003 
[5] . Corollary model dependent quests for a candidate particle have been carried 
out in some of the many possible model frameworks and have been improved 
with time. In particular, some scenarios either for mixed spin- independent (SI) 
and spin-dependent (SD) coupled WIMPs or for purely SI coupled WIMPs or 
for purely SD ^ coupled WIMPs have been considered in some of the many 
possible model frameworks as well as the case of WIMPs with preferred inelastic 
scattering. 

Finally, some comments on the claims for contradiction by Cdms-I, Edel- 
weiss-I, Zeplin I can be found in ref. [5]^, where some recent positive hints - 
not in contradiction with the DAMA/Nal result - from Dark Matter indirect 
searches are also summarised. 

3 Final DAMA/Nal model independent result on WIMP 
annual modulation signature over 7 annual cycles 

A model independent investigation of the WIMP annual modulation signature 
has been realised by exploiting the time behaviour of the experimental residual 

rates of the single-hit events in the lowest energy region over seven annual cy- 
cles (total exposure: 107731 kg x day) [5], as performed in the past on partial 
exposures. 

The exposure of each annual cycle, DAMA/NaI-1 to -7, is reported in Table 
1, where the exposure of the DAMA/NaI-0 running period (partially overlapped 
with DAMA/Nal- 1) is quoted as well. 

This model independent approach on the data of the seven annual cyc;lcs 
offers an immediate evidence of the presence of an annual modulation of the 
rate of the single-hit events in the lowest energy region as shown in Fig. 1 - left, 
where the time behaviours of the measured (2-4), (2-5) and (2-6) keV single-hit 
events residual rates are depicted. 

^ We remind that JHEP 0107 (2001) 044 is not at all in conflict with a purely SD 
solution since it considered only two particular purely SD couplings (of the many 
possible) in a strongly model dependent context and using modulation amplitudes 
valid instead only in a particular purely SI case. Moreover, the mixed SI & SD case 
was not involved at all in that discussion. 

* Here we just remind e.g. that no model independent comparison is possible among 
those experiments and DAMA/Nal because of the different methodological ap- 
proaches, of the different target nuclei, etc.. As regards possible model dependent 
comparisons, those experiments give result in a single purely SI model framework 
with fixed/selected assumptions, neglecting experimental and theoretical uncertain- 
ties and - as also done in this conference - have quoted so far the DAMA/Nal result 
in an incorrect , unupdated and incomplete way and they have been and axe ignoring 
the existence of other solutions. 
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Table 1. Exposure of each annual cycle (DAMA/Nal- 1 to -7) and of the DAMA/NaI-0 
running period [6,17-24,5]. Sec text. 

I i-'cri.,(ls I Exposure (kg X day) | 



|DAMA/Nal-0| 4123.2 



DAMA/NaI-1 


3363.8 (winter) 
+ 1185.2 (summer) 


DAMA/NaI-2 


14962 
(Nov. — > end of July) 


DAMA/NaI-3 


22455 

(middle Aug. — ^ end of Sept.) 


DAMA/NaI-4 


16020 

(middle Oct. — middle Aug.) 


DAMA/NaI-5 


15911 
(Aug. — ^ end of July) 


DAMA/NaI-6 


16608 
(Nov. end of July) 


DAMA/NaI-7 


17226 
(Aug. — > end of July) 


Total 1-7 


107731 



The data favour the presence of a modulated cosine-hke behaviour {A- cosa;(t— 
to)) at 6.3 a C.L. and their fit for the (2-6) kcV larger statistics energy interval of- 
fers modulation amplitude equal to (0.0200±0.0032) cpd/kg/keV, to = (140±22) 
days and T = ^ = (1.00 ± 0.01) year, all parameters kept free in the fit. The 
period and phase agree with those expected in the case of a WIMP induced effect 
(T = 1 year and to roughly at ~ 152.5*'' day of the year). The test on the 
(2-6) keV residual rate in Fig. 1 - left disfavours the hypothesis of unmodulated 
behaviour giving a probability of 7 • 10~^ {x^/d.o.f. = 71/37). We note that, 
for simplicity, in Fig. 1 - left the same time binning already considered in ref. 
[20,21] has been used; the result of this approach is similar by choosing other 
time binnings. The experimental residuals given in Fig. 1 - left have also been 
fitted, according to the previous procedure, fixing the period at 1 year and the 
phase at 2""^ June; the best fitted modulation amplitudes are: (0.0233 ± 0.0047) 
cpd/kg/keV for the (2-4) keV energy interval, (0.0210 ±0.0038) cpd/kg/keV for 
the (2-5) keV energy interval, (0.0192 ± 0.0031) cpd/kg/keV for the (2-6) keV 
energy interval, respectively. 

The same data have also been investigated by a Fourier analysis (performed 
according to ref. [25] including also the treatment of the experimental errors and 
of the time binning), obtaining the result shown in Fig. 1 - right, where a clear 
peak corresponding to a period of 1 year is evident. 

In Fig. 2 the experimental single-hit residual rate from the total exposure 
of 107731 kg X day is presented, as in a single annual cycle, for two different 
energy intervals; as it can be seen the modulation is clearly present in the (2-6) 
keV energy region, while it is absent just above. 

Finally, a suitable statistical analysis has shown that the modulation ampli- 
tudes are statistically well distributed in all the crystals, in all the data taking 
periods and considered energy bins [5]. 

On the investigation of possible systematic effects and side reactions. 

As previously mentioned, to mimic the annual modulation signature a system- 
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Fig. 1. On the left: experimental residual rate for single-hit events in the (2-4), (2-5) 
and (2-6) keV energy intervals as a function of the time over 7 annual cycles (total 
exposure 107731 kg x day); end of data taking July 2002. The experimental points 
present the errors as vertical bars and the associated time bin width as horizontal bars. 
The superimposed curves represent the cosinusoidal functions behaviours expected for 
a WIMP signal with a period equal to 1 year and phase at 2"'' June; the modulation 
amplitudes have been obtained by best fit [5]. On the right: power spectrum of the 
measured (2-6) keV single-hit residuals calculated according to rcf. [25] , including also 
the treatment of the experimental errors and of the time binning. As it can be seen, 
the principal mode corresponds to a frequency of 2.737 • 10~^ d~^, that is to a period 
of ~ 1 year. 



atic effect or side reaction should not only be able to account for the whole 
observed modulation amplitude, but also able to satisfy the requirements of a 
WIMP induced effect. A careful investigation of all the known possible sources 
of systematics and side reactions has been regularly carried out and published 
at time of each data release. In particular, detailed quantitative discussions can 
be found in ref. [21,5] and will not be repeated here ^. As it can be seen there, 



^ We take this opportunity only to comment that the sizeable discussions reported e.g. 
in [21,5] already demonstrated that a possible modulation of neutron flux (possibly 
observed by the ICARUS coll. as reported in the ICARUS internal report TM03-01) 
cannot quantitatively contribute to the DAMA/Nal observed modulation amplitude, 
even if the neutron flux would be assumed to be 100 times larger than measured 
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Fig. 2. Experimental single-hit residual rate from the total exposure of 107731 kg x day 

as in a single annual cycle. The experimental points present the errors as vertical bars 

and the associated time bin width as horizontal bars. The initial time is taken at August 

7*''. Fitting the data with a cosinusoidal function with period of 1 year and phase at 

152.5 days, the following amplitudes are obtained: (0.0195 ± 0.0031) cpd/kg/keV and 

— (0.0009 ± 0.0019) cpd/kg/keV, respectively. Thus, a clear modulation is present in 

the lowest energy region, while it is absent just above. 



no systematic effect or side reaction able to mimic a WIMP induced effect has 
been found. 

A further result: the multiple-hits events. As a further relevant investigation, 
the multiple-hits events also collected during the DAMA/NaI-6 and 7 running 
periods (when each detector was equipped with its own Transient Digitizer with 
a dedicated renewed electronics) have been studied and analysed by using the 
same identical hardware and the same identical software procedures as for the 
case of the single-hit events. The multiple-hits events class - on the contrary 
of the single-hit one does not include events induced by WIMPs since the 
probability that a WIMP scatters off more than one detector is negligible. 

Fig. 3 shows the behaviour of the residual rate of multiple-hits events in the 
(2-6) keV energy interval measured during the DAMA/NaI-6 and -7 running 
periods as a function of the time in a year. It is compared with the residual rate 
of the single-hit events measured in the same energy interval with the total expo- 
sure. Fitting these data with the function A- cosLv{t—to) with period of 1 year and 
phase at 152.5 days, the following amplitudes are obtained: A = (0.0195±0.0031) 
cpd/kg/keV and A = —(3.9 ± 7.9) • 10~^ cpd/kg/keV for single-hit and multiple- 
hits residual rates, respectively. Thus, a 6.3 a C.L. evidence of annual modulation 
is present in the single-hit residuals (events class to which the WIMP-induced 
recoils belong), while it is absent in the multiple-hits residual rate (event class to 
which only background events belong). Since the same identical hardware and 
the same identical software procedures have been used to analyse the two classes 
of events, the obtained result offers an additional strong support for the presence 
of Dark Matter particles in the galactic halo further excluding any side effect 
either from hardware or from software procedures or from background. 



at LNGS by several authors with different techniques over more than 15 years; in 
addition, as widely known, it cannot satisfy all the peculiarities of the signature. 
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Fig. 3. Experimental residual rates over seven annual cycles for single-hit events (open 
circles) - class of events to which WIMP events belong - and over the last two annual 

cycles for multiple-hits events (filled triangles) class of events to which WIMP events 
do not belong - in the (2-6) keV cumulative energy interval. They have been obtained 
by considering for each class of events the the data as collected in a single annual cycle 
and using in both cases the same identical hardware and the same identical software 
procedures. The initial time is taken on August 7*^. 

Conclusion. In conclusion, the presence of an annual modulation in the single- 
hit events residual rate in tlie lowest energy interval (2 - 6) keV satisfying all the 
features expected for a WIMP component in the galactic halo is supported by 
the data of the seven annual cycles at 6.3 a C.L.. This is the experimental result 
of DAMA/Nal. It is model independent; no other experiment whose result can 
be directly compared with this one is available so far in the field of Dark Matter 
investigation. 



4 Corollary result: quests for a candidate peirticle in 
some model frameworks with the data of the seven 
annual cycles 

On the basis of the obtained model independent result, corollary investigations 
can also be pursued on the nature and coupling of the WIMP candidate. This 
latter investigation is instead model dependent and - considering the large un- 
certainties which exist on the astrophysical, nuclear and particle physics assump- 
tions and on the parameters needed in the calculations - has no general meaning 
(as it is also the case - as mentioned above - of exclusion plots and of the WIMP 
parameters evaluated in indirect search experiments). Thus, it should be han- 
dled in the most general way as we have preliminarily pointed out with time in 
the past [17-24] and we have discussed in some specific details in ref. [5]; other 
efforts on this topic are also in progress. Candidates, kinds of WIMP couplings 
with ordinary matter and implications, cross sections, nuclear form factors, spin 
factors, scaling laws, halo models, priors, etc. are discussed in details in ref. [5] 
and we invite the reader to this reference since these arguments arc necessary to 
correctly understand the results obtained in corollary quests and the real validity 
of any claimed model dependent comparison in the field. The results presented in 
ref. [5] and summarised here are, of course, not exhaustive of the many possible 
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scenarios which at present level of knowledge cannot be disentangled. Some of 

the open questions are: i) which is the right nature for the WIMP particle ^; 
ii) which is its right couplings with ordinary matter (mixed SI&SD, purely SI, 
purely SD or preferred inelastic) iii) which are the right form factors and related 
parameters for each target nucleus; iv) which is the right spin factor for each 
target nucleus (some nuclei are disfavoured to some kinds of interactions; for ex- 
ample, in case of an interaction with SD component even a nucleus sensitive in 
principle to SD interaction could be blinded by the spin factor if unfavoured by 
the 9 value^); v) which are the right scaling laws (let us consider as an example 
that even in a MSSM framework with purely SI interaction the scenario could be 
drastically modified as discussed recently in ref. [27]); vi) which is the right halo 
model and related parameters; vii) which are the right values of the experimental 
parameters within their uncertainties; etc. As an example, we remind that not 
only large differences in the measured rate can be expected when using target 
nuclei sensitive to the SD component of the interaction (such as e.g. ^^iVa and 
^^^/) with respect to those largely insensitive to such a coupling (such as e.g. 
""■*Ge and "°*S'i), but also when using different target nuclei although all - in 
principle - sensitive to such a coupling (compare e.g. the Xenon and Tellurium 
cases with the Sodium and Iodine cases) [5]. 

In the following some of the results discussed for some of the many possi- 
ble model dependent quests for a WIMP candidate are briefly reminded [5]. In 
particular, they have been obtained from the data collected during all the seven 
annual cycles, considering the halo models summarized in [24,5] for three of the 
possible values of the local velocity Vq: 170 km/s, 220 km/s and 270 km/s. The 
escape velocity has been maintained at the fixed value: 650 km/s. It is worth to 
note that the present existing uncertainties on the knowledge of the escape ve- 
locity can play a relevant role in evaluating allowed regions (and corresponding 
best fit values for WIMP mass and cross section) e.g. in the cases of preferred in- 



Several candidates fulfil the cosmological and particle Physics requirements neces- 
sary in order to be considered as a Dark Matter candidate particle: not only the 

neutralino foreseen in the suporsymmctric theories, but also a heavy neutrino of a 
4"* family (there is still room for it as reported in literature), the sneutrino in the 
scenario described by [14] (providing - through the transition from lower to up- 
per mass eigenstate - preferred inelastic scattering with target-nuclei), the "mirror" 
Dark Matter [26], etc.. Moreover, in principle whatever Weakly Interacting, neutral, 
(quasi-)stable and Massive (whose acronym is WIMP) particle, even not yet foreseen 
by a theory, can be a suitable candidate. As regards in particular the neutralino, we 
note that the theories have not stringent predictive capability for its cross sections 
and for its mass because of the large number of free parameters in the theory and 
of the several assumptions required; thus, e.g. the expectations for its nuclear cross 
sections span over several orders of magnitude as it can also be seen in literature. In 
addition, we take this occasion to remind that the neutralino has both SI and SD 
couplings with the ordinary matter. 

We remind that tgO = an/ap is the ratio between the WIMP-neutron and the WIMP- 
proton effective SD coupling strengths. On and Op, respectively [22,5]; 6 is defined in 
the [0,7r) interval. 
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elastic WIMPs and of light mass WIMP candidates; its effect would be instead 
marginal at large WIMP masses. The possible scenarios have been exploited 
for those halo models in some discrete cases including some of the uncertain- 
ties which exist in the parameters of the used nuclear form factors and in the 
quenching factors; for the details see ref. [5]. The results summarised here are 
not exhaustive of the many scenarios possible at present level of knowledge: e.g. 
for some other recent ideas see the already quoted [15,27]. 

For simplicity, here the results of these corollary quests for a candidate par- 
ticle are presented in terms of allowed regions obtained as superposition of the 
configurations corresponding to likelihood function values distant more than 4f7 
from the null hypothesis (absence of modulation) in each of the several (but still 
a limited number) of the possible model frameworks considered here. Obviously, 
larger sensitivities than those reported in the following figures would be reached 
when including the effect of other existing uncertainties on the astrophysical, 
nuclear and particle Physics assumptions and related parameters; similarly, the 
set of the best fit values would also be enlarged as well. 

As well known, DAMA/Nal is intrinsically sensitive both to low and high 
WIMP mass having both a light (the ^^Na) and a heavy (the ^^'^I) target- 
nucleus; in previous corollary quests WIMP masses above 30 GeV (25 GeV in 
ref. [17]) have been presented [18,20,22-24] for few (of the many possible) model 
frameworks. However, that bound holds only for neutralino when supersymmet- 
ric schemes based on GUT assumptions are adopted to analyse; the LEP data 
[28]. Thus, since other candidates are possible and also other scenarios can be 
considered for the neutralino itself as recently pointed out the present model 
dependent lower bound quoted by LEP for the neutralino in the supersymmetric 
schemes based on GUT assumptions (37 GeV [32]) is simply marked in the fol- 
lowing figures. It is worth to note that when this mass limit is adopted, it selects 
the WIMP-Iodine elastic scattering as dominant because of the used scaling laws 
and of kinematical arguments. Finally, the prior from DAMA/NaI-0 has properly 
been considered as well. 

WIMPs with mixed SI&SD interaction. The most general scenario of 
WIMP nucleus elastic interaction, to which the DAMA/Nal target nuclei are 
fully sensitive, is the one where both the SI and the SD components of the cross 
section are present. Thus, as first we introduce here the case for a candidate 
with both SI and SD couplings to ordinary matter [5] similarly as we did in the 
past also in ref. [22] on partial exposure. In this general scenario the space of the 
free parameters is a 4-dimensional volume defined by m.w, S,<^si ^<^sd and 
(which varies from to tt). Thus, the general solution would be a four dimen- 
sional allowed volume for each considered model framework. Since the graphic 
representation of this allowed volume is quite difficult, we only show in Fig. 4 
the obtained regions in the plane £,asi vs ^aso for some of the possible and 



® In fact, wlien the assumption on tlie gaugino-mass unification at GUT scale is re- 
leased neutralino masses down to ~ 6 GeV are allowed [29-31]. 
® ^ (^ < 1) is defined here as the fractional amount of local WIMP density. 
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Fig. 4. Case of a WIMP with mixed SI&SD interaction in the given model frameworks. 
Coloured areas: example of slices (of the allo-wed volume) in the plane £,0-31 vs £,asD 
for some of the possible mw and 9 values. Inclusion of other existing uncertainties on 
parameters and models would further extend the regions; for example, the use of more 
favourable form factors and/or of more favourable spin factors than the considered 
ones -would move them towards lower cross sections. For details see [5]. 



mw values in the model frame"works considered here. In particular, we report 
just four couplings, -which correspond to the follo-wing values of the mixing angle 
9: i) 9 = (a„ =0 and Op 7^ or \ap\ >> |a„|) corresponding to a particle -with 
null SD coupling to neutron; ii) 9 = tt/A (up = Un) corresponding to a particle 
"with the same SD coupling to neutron and proton; iii) 9 = Tr/2 (a„ ^ and ap = 
or |a„| >> I Op I) corresponding to a particle -with null SD couplings to proton; 
iv) 9 = 2.435 rad — -0.85) corresponding to a particle -with SD coupling 

through Zq exchange. The case ap = — a„ is nearly similar to the case iv). 

From the given figures it is clear that at present either a purely SI or a purely 
SD or a mixed SI&SD configurations are supported by the experimental data of 
the seven annual cycles. 

WIMPs -with dominant SI interaction. Generally, the case of purely SI cou- 
pled WIMP is mainly considered in literature. In fact, often the spin-independent 
interaction -with ordinary matter is assumed to be dominant since e.g. most of 
the used target-nuclei are practically not sensitive to SD interactions (as on the 
contrary ^^Na and ^^^I are) and the theoretical calculations are even much more 
complex and uncertain. 
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Thus, following an analogous procedure as for the previous case, we have 
exploited for the same model frameworks the purely SI scenario. In this case the 
free parameters are two: mw and ^asi- 
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Fig. 5. On the left : Case of a WIMP with dominant SI interaction for the given 
model frameworks. Region allowed in the plane {mw, (,<^si)- The vertical dotted line 
represents a bound in case of a neutralino candidate when supersymmetric schemes 
based on GUT assumptions are adopted to analyse the LEP data; the low mass region 
is allowed for neutralino when other schemes are considered and for every other WIMP 
candidate; see text. While the area at WIMP masses above 200 GeV is allowed only 
for few configurations, the lower one is allowed by most configurations (the colored 
region gathers only those above the vertical line) [5]. The inclusion of other existing 
uncertainties on parameters and models would further extend the region; for example, 
the use of more favourable SI form factor for Iodine alone would move it towards lower 
cross sections. On the right : Example of the effect induced by the inclusion of a SD 
component different from zero on allowed regions given in the plane ^crsi vs mw- In 
this example the Evans' logarithmic axisymmetric C2 halo model with vq = 170 km/s, 
po equal to the maximum value for this model and a given set of the parameters' values 
(see [5]) have been considered. The different regions refer to different SD contributions 
for the particular case of = 0: asD = pb (a), 0.02 pb (b), 0.04 pb (c), 0.05 pb (d), 
0.06 pb (e), 0.08 pb (f). Analogous situation is found for the other model frameworks. 



In Fig. 5 - left the region allowed in the plane mw and ^asi for the consid- 
ered model frameworks is reported. The configurations below the vertical line 
are of interest for neutralino when the assumption on the gaugino-mass unifica- 
tion at GUT scale is released and for every other kind of WIMP candidate. As 
shown in Fig. 5 - left, also WIMP masses above 200 GeV are allowed for some 
configurations; details can be found in ref. [5]. Of course, best fit values of cross 
section and WIMP mass span over a large range in the considered frameworks. 

Let us now point out, in addition, that configurations with ^crsi even much 
lower than those shown in Fig. 5 - left are accessible in case an even small SD 
contribution is present in the interaction. This possibility is clearly pointed out 
in Fig. 5 - right where an example of allowed regions in the plane {mw, ^o'si) 
corresponding to different SD contributions is reported for the case 6* = 0. As it 
can be seen, increasing the SD contribution the allowed regions involve SI cross 
sections much lower than 10~^ pb. It can be noted that for asD > 0.08 pb the 
annual modulation effect observed is also compatible - for mw — 40 — 75 GeV 
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- with a WIMP candidate with no SI interaction at alL Analogous situation is 
found for the other model frameworks. 

WIMPs with dominant SD interaction. Let us now focus on the case of 
a candidate with purely SD coupling to which DAMA/Nal is - as mentioned - 
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Fig. 6. Case of a WIMP with dominant SD interaction in the given model frameworks. 
An example of the region allowed in the plane (mw, ^crso); here 6 = 2.435, Zo coupling, 
(9 is defined in the [0, tt) range). For the definition of the vertical line and of the 
coloured area see previous figure caption; see also text. Inclusion of other existing 
uncertainties on parameters and models (as discussed in ref. [5]) would further extend 
the SD allowed regions. For example, the use of more favourable SD form factors and/or 
more favourable spin factors would move them towards lower cross sections. Values of 
^(TSD lower than those corresponding to these allowed regions are possible also e.g. in 
case of an even small SI contribution (see ref. [5]). 

fully sensitive. 

When the SD component is different from zero, a very large number of possi- 
ble configurations is available. In fact, in this scenario the space of free parame- 
ters is a 3-dimensional volume defined by mw, £,<Jsd and 6 (which can vary from 
to tt). Here, just as an example we show the results obtained only for a partic- 
ular coupling, which correspond to a mixing angle 9 = 2.435 {Zq coupling); see 
Fig. 6; other configurations are possible varying the 9 value. The area at WIMP 
masses above 200 GeV is allowed for low local velocity and all considered sets of 
parameters by the Evans' logarithmic C2 co- rotating halo model [5]. Moreover, 
the accounting for the uncertainties e.g. on the form factors and/or on the spin 
factors as well as different possible formulations of the SD form factors would 
extend the allowed regions, e.g. towards lower (,asD values. Finally, £,asD lower 
than those corresponding to the regions shown in Fig. 6 are possible also e.g. in 
case of an even small SI contribution (see ref. [5]). 

WIMPs with preferred inelastic interaction. An analysis considering the 
same model frameworks has been carried out for the case of WIMPs with pre- 
ferred inelastic interaction [14] as we did also in the past in ref. [23] on partial 
exposure. 

In this inelastic Dark Matter scenario an allowed volume in the space (fcTp, 
mw, 5) is obtained; 5 is the mass splitting of the WIMP particle which can be 
excited following an inelastic interaction [14,23,5]. For simplicity, Fig. 7 shows 
slices of such an allowed volume at some given WIMP masses. There the super- 
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Fig. 7. Case of a WIMP with preferred inelastic interaction in the given model frame- 
works. Examples of slices (coloured areas) of the allowed volumes (^fp, 5, raw) for 
some mw values for WIMP with preferred inelastic interaction. Inclusion of other ex- 
isting uncertainties on parameters and models would further extend the regions; for 
example, the use of a more favourable SI form factor for Iodine and different escape 
velocity would move them towards lower cross sections [5]. 

positions of the allowed regions obtained, when varying the model framework 
within the considered set [5], are shown for each tow We remind that in these 
calculations Vesc has been assumed at fixed value (as in the previous cases) , while 
its present uncertainties can play a significant role in this scenario of WIMP with 
preferred inelastic scattering. 

Conclusion. In this section we have briefly summarized some quests for the 
candidate particle in some of the many possible scenarios. We further stress 
that, although several scenarios have been investigated, these corollary analyses 
are not exhaustive at all because of the present poor knowledge on many astro- 
physical, nuclear and particle Physics needed assumptions; moreover, additional 
scenarios can also be possible as also shown e.g. by some recent papers appeared 
in literature. Other model dependent quests are already under consideration. 

5 The second generation: DAMA/LIBRA 

In 1996 a ton set-up was proposed by Bernabei et al.; as a consequence a new 
R&D for Nal(Tl) radiopurification has been carried out and the second genera- 
tion set-up DAMA/LIBRA (~250 kg Nal(Tl)) has been funded and reahsed as 
an intermediate step. This R&D with Crismatec-St. Gobain company has ex- 
ploited for the first time new chemical/physical radiopurification procedures in 
Nal and Til powders already selected for radiopurity. In addition, new selected 
materials and set-up components as well as new protocols have been employed 
for building DAMA/LIBRA. 
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In 2002 ended the production of detectors and of new parts of the mstaUation 
for DAMA/LIBRA. Thus, after July 2002 - at the completion of its data taking - 
DAMA/Nal was fully dismounted and the installation of the new DAMA/LIBRA 
started. The experimental site as well as many components of the installation 
itself were implemented (environment, shield of PMTs, wiring, HP Nitrogen 
system, cooling water of air conditioner, electronics and DAQ, etc.). Before 
the installation, all the Cu parts were chemically etched following a devoted 
protocol and maintained in HP Nitrogen atmosphere. In addition, all the proce- 
dures performed during the dismounting of DAMA/Nal and the installation of 
DAMA/LIBRA were carried out in HP Nitrogen atmosphere. This was realised 
by using a Scuba system (a self-contained underwater breathing apparatus) mod- 
ified in order to avoid that the entire breath is expelled into the surrounding air 
when the operator exhales; the air cylinders were kept five meters away and the 
output line was two meters long. 



Fig. 8. DAMA/LIBRA: a view of the detectors inside the inner Cu box at end of 
installation (photo taken in HP Nitrogen atmosphere). All the used materials have 
been deeply selected for radiopurity. 

The new DAMA/LIBRA, having an exposed mass of ~ 250 kg, an higher 
overall radiopurity and improved performances, offers an increased experimental 
sensitivity to further investigate the DAMA/Nal observed effect and to improve 
investigations on the nature of the candidate particle trying to disentangle among 
different possible astrophysical, nuclear and particle Physics models as well as 
other new possible scenarios. See just as an example of new ideas the case of the 
mirror Dark Matter [26] , the case of a contribution to the dark halo from the 
Sagittarius Dwarf Tidal Stream [15], the case of possible different nuclear scaling 
laws even for the neutralino candidate in MSSM [27], etc.. Moreover, the low 
background DAMA/LIBRA offers a powerful tool for the Dark Matter investi- 
gation in the future since it is e.g. sensitive: i) both to low (through interaction 
on ^^Na) and to high (through interaction on ^^^I) mass Dark Matter particles; 
ii) both to mixed SI&SD, to purely SI, to purely SD couplings and to preferred 
inelastic scattering as well as to other possible kind of Dark Matter candidates 
(e.g. mirror Dark Matter). 

DAMA/LIBRA has started the preliminary data taking on March 2003 and 
it has been planned to run for several years. At present a new R&D effort toward 
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the possible Nal(Tl) ton set-up has been funded and related works have already 
been started. 

6 Conclusion 

DAMA/Nal has been a pioneer experiment running at LNGS for about a decade 
and investigating as first the WIMP annual modulation signature with suitable 
exposed mass, sensitivity and control of the running parameters. During seven 
independent experiments of one year each one, it has pointed out the presence 
of a modulation satisfying the many peculiarities of a WIMP induced effect, 
reaching an evidence at 6.3 cr C.L.. No other experiment has so far been in 
position to give a result directly comparable in a model independent way with 
that of DAMA/Nal. As a corollary, it has also pointed out the complexity of 
the quest for a WIMP candidate because of the present poor knowledge on the 
many astrophysical, nuclear and particle physics aspects. 

After the completion, on July 2002, of the DAMA/Nal data taking, the sec- 
ond generation DAMA/LIBRA set-up (~ 250 kg mass) has been installed and 
preliminarily put in operation, as a result of continuous efforts by the DAMA col- 
laboration toward the creation of ultimate radiopure Nal(Tl). DAMA/LIBRA, 
having a larger exposed mass and an higher overall radiopurity, will signifi- 
cantly contribute in the incoming years to the further understanding of the 
field. Presently, new R&D efforts toward a possible Nal(Tl) ton set-up have 
been funded and related works have already been started. 
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